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' We investigate extremal dyon black holes in the Einstein-Maxwell-dilaton (EMD) theory with 

' higher curvature corrections in the form of the Gauss-Bonnet density coupled to the dilaton. In 

the same theory without the Gauss-Bonnet term the extremal dyon solutions exist only for discrete 
values of the dilaton coupling constant a. We show that the Gauss-Bonnet term acts as a dyon 
rzi ' hair tonic enlarging the allowed values of a to continuous domains in the plane (a, qm) the second 

parameter being the magnetic charge. In the limit of the vanishing curvature coupling (a large 
magnetic charge) the dyon solutions obtained tend to the Reissner-Nordstrom solution but not to 
the extremal dyons of the EMD theory. Both solutions have the same values of the horizon radius as 
a function of charges. The entropy of new dyonic black holes interpolates between the Bekenstein- 
Hawking value in the limit of the large magnetic charge (equivalent to the vanishing Gauss-Bonnet 
coupling) and twice this value for the vanishing magnetic charge. Although an expression for the 
entropy can be obtained analytically using purely local near-horizon solutions, its interpretation as 
the black hole entropy is legitimate only once the global black hole solution is known to exist, and 
we obtain numerically the corresponding conditions on the parameters. Thus, a purely local analysis 
is insufficient to fully understand the entropy of the curvature corrected black holes. We also find 
dyon solutions which are not asymptotically flat, but approach the linear dilaton background at 
infinity. They describe magnetic black holes on the electric linear dilaton background. 

> 

O ' PACS numbers: 04.20.Jb, 04.65. +e, 98.80.-k 

(N ■ 



I. INTRODUCTION 



0^ 
O 
00 . 

' In a recent paper [l| we started an investigation of extremal black holes in the four-dimensional Einstein-Maxwell- 
dilaton-Gauss-Bonnet model motivated by an interest to microscopic string calculations of the black hole entropy (for a 
review see @, Q)- In the theories with higher curvature corrections, the entropy deviates from the Bekenstein-Hawking 
/\ ' value and can be calculated using Wald's formalism 0, H, H, 0, B @ ■ Remarkably, it still exhibits an agreement with the 
string theory predictions at the corresponding level, both in the BPS [H [11 B [11, [11, [IS H [11 S 11^^ 



- - ' and non-BPS [H, [H [H, [H, [H [H, [^, [H [Sll, HI, [H [M [H, [H, HSlil,!!!! cases. In some supersymmetric models 
with higher curvature terms exact classical solutions for static black holes were obtained [l8, 19, 22]. Moreover, as 
was argued by Sen [iO, [HI [13 , the knowledge of the global black hole solutions is not necessary in order to compare 
classical and quantum results for the entropy: in the classical theory the entropy can be computed locally using the 
entropy function approach based on the attractor property typical for supergravity black holes [11[23, [25, 26, 27, 28,]. 
However the question remains, whether the global black holes corresponding to local solutions used to construct the 
entropy function do really exist. Generically, the existence of local solutions exhibiting the event horizons does not 
guarantee that they describe black holes which must be regular outside the horizon and asymptotically flat. Even 
though formally the entropy can be obtained from local considerations, the parameters involved may be subject to 
restrictions which are revealed only when we try to extend the local solutions to infinity. This issue is addressed in 
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the present paper for extremal dyon black holes with the AdS2 x S'^ horizon in the Gauss-Bonnet gravity. 

The extremal dilatonic black hole [i^ 0, is a particularly interesting model associated with the heterotic 

string. In the Einstein-Maxwell-dilaton (EMD) theory without curvature corrections such a solution have a singular 
horizon of zero radius, the corresponding Bekenstein-Hawking entropy being zero. Typical higher curvature correction 
to this theory is given by the Gauss-Bonnet density coupled to the dilaton (later on referred as the Einstein-dilaton- 
Gauss-Bonnet (EDGB) model). In this model the local solutions with the AdS2 x (Reissner-Nordstrom type) 
horizons of finite radius can be found. Presumably they should describe the extremal black holes possessing non-zero 
entropy. In Q we have explored whether the local solutions constructed as series expansions in the vicinity of the event 
horizon can be extended to infinity as asymptotically flat black holes. It was shown that this is possible only within 
some bounded region in the space of parameters. Actually, the model contains two parameters which are worth to be 
considered as independent in the classical theory: the Gauss-Bonnet coupling constant a weighting the contribution 
of the Gauss-Bonnet term, and the dilaton coupling constant a. In the case of the purely electric configurations, the 
local solutions around the AdS2 x horizon exist only for non-zero a, since no dilatonic purely electric (or purely 
magnetic) black holes with the AdS2 x S*^ horizon are possible in the Einstein-Maxwell-dilaton (EMD) theory without 
curvature corrections. Thus the only continuously varying parameter is the dilaton coupling constant a. If a = 0, the 
model reduces to the Einstein-Maxwell theory, in which the desired extremal solution do exist (the extremal Reissner- 
Nordstrom solution). We have studied whether similar solutions exist in the curvature-corrected theory for non-zero 
a, aiming to investigate the case a — 1 relevant to the compactified heterotic string. It turned out that the extremal 
black hole solutions exist only in a bounded region < a < where the critical value is close to 1/2. For greater 
a, the local solutions exhibiting the AdS2 x horizons can not be extended to infinity as regular black holes, but 
develop singularities outside the horizon. The threshold value of a might serve an indication that the corresponding 
string configuration experiences some qualitative change of state like the black hole-string transition (47j . 

The case of the purely electric extremal black hole in the EDGB model is special in the sense that no similar 
solutions exist in the EMD theory without curvature corrections (unless a = 0), so there is no smooth transition to 
the non-corrected theory. The situation becomes more flexible if we allow for both electric and magnetic charges to be 
present, since extremal dyons with the AdS2 x horizon do exist in the EMD theory. As we will see, such solutions 
are possible for discrete values of the dilaton coupling constant: = 1 (the heterotic case), — 3 (the Kaluza-Klcin 
case) and some sequence of other integer (coinciding with the sequence found in [gJ] for existence of non-extremal 
black holes with two horizons). Therefore in the curvature-corrected theory we have an infinite sequence of starting 
points for the dilaton coupling constant, not just the trivial point a = 0. So one can expect to have much larger 
domain of existence of extremal charged dyonic black holes in the EDGB model than in the purely electric case. 

We will be interested here uniquely by the extremal black holes with the degenerate event horizon. Non-degenerate 
black holes in the same theory were extensively studied in the past both perturbatively [H, and numerically [sol . 
[sTI. [5^ [53l [53 . More recently global properties of EDGB black hole solutions were studied using the dynamical system 
approach [53, HE 113, Hll . Stability issues were discussed in [H, [13, HH, [sl, [G^I . In these papers the existence of both 
neutral and charged asymptotically flat solutions with a non-degenerate event horizon and without naked singularities 
was established. These solutions have the Schwarzschild type event horizon and they do not possess the extremal 
limits. The solution with the degenerate event horizon thus form a separate branch of EDGB black holes which was 
not studied before. 

The existence of non-degenerate dyonic black holes in the EMD theory without curvature corrections was studied in 
detail in J34j. It was found that such solutions generically exhibit one (non-degenerate) horizon of the Schwarzschild 
type, but for some discrete values of the dilaton coupling constant there are solutions with two horizons of the Reissner- 
Nordstrom type. In this latter case the limiting extremal solutions turn out to be possible. We will show that adding 
the Gauss-Bonnet term to the EMD action acts as a hair tonic for extremal dyons, allowing for continuously varying 
dilaton couplings. Still, the allowed domain is bounded (in somewhat irregular way), so the threshold behavior 
observed in pj] persists in the dyon case too. The solutions corresponding to the values of parameters approaching 
the boundary of the allowed domain exhibit interesting saturation properties similar to the BPS bounds in the EMD 
theory without curvature corrections. 

The plan of the paper is as follows. The Sec. 2 contains general definitions and investigation of the hidden symmetries 
of the reduced theory. It is shown that the one-dimensional theory possesses the two-parametric off-shell and the 
three-parametric on-shell symmetry groups. These symmetries serve as a convenient tool allowing to describe the 
space of solution semi-analytically, in spite of the fact that the equations of motion do not have analytic solutions. In 
Sec. 3 we investigate local series solutions around the assumed AdS2 x S*^ horizon. We show that for fixed values of 
the electric and magnetic charges the local solutions contain only one free parameter, contrary to two parameters in 
the corresponding theory without curvature corrections. Examining the higher order coefhcients we find the family 
of curves in the parameter plane of a and qm (the magnetic charge) on which the coefficients of the local solution 
are singular. These curves are parameterized by the same sequence of integers af which correspond to existence of 
dyons in the uncorrected theory. In this section we also calculate the entropy and show that it interpolates between 
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the Bekenstein-Hawking value A/ 4 in the case of vanishing Gauss-Bonnet couphng (or large magnetic charge) and 
twice the Bekenstein-Hawking large value A/2 for purely electric solutions, as found in The Sec. 4 is devoted to 
asymptotic solutions. We find that, similarly to the uncorrected EMD theory, there are two physically interesting 
asymptotic patterns corresponding either to usual asymptotically flat black holes, or to black holes on the linear 
dilaton background (65j . Numerical results are presented in the Sec. 5. Using the on shell symmetries we can express 
the ratios of physical parameters as functions of the dilaton coupling constant only, which clarifies the properties of the 
whole family of dyonic solutions obtained. We also demonstrate their asymptotic BPS-type behavior on the boundary 
of the allowed domain of parameters. The parameter region for the asymptotically LDB solutions is shown to locate 
on the other side of the limiting singular curve in parameter plane. In the Appendix A we give some details concerning 
the EMD dyons and show how analytic solutions known for two lower values of the discrete dilaton coupling sequence 
can be obtained by summing up the local series solutions valid in the vicinity of the horizon. In the Appendix B 
some higher order coefficient of the local solution in the EDGB model are listed which are necessary in deriving the 
existence of limiting curves in the parameter space. 



II. GENERAL SETTING 



We consider the four-dimensional dilatonic Gauss-Bonnet theory (EDGB) which is the Einstein-Maxwell-dilaton 
theory (EMD) with an arbitrary dilaton coupling constant a modified by the Gauss-Bonnet (GB) term: 

^ = / - '^Md"^ - e'"* {F^ - cCgb) } V^d^x, (1) 
where £gb is the Gauss-Bonnet density 

/:GB - i?' - 4i?^.i?^'^ + i?a/3M-^i?"''^^ (2) 

This action contains two parameters (we use the units G = c = 1): the dilaton coupling constant a and the Gauss- 
Bonnet coupling constant a. We assume a > 0, a > 0, solutions for negative a can be obtained changing the sign 
of the dilaton. Note that in this action the Maxwell term is not multiplied by a to facilitate decoupling of the 
Gauss-Bonnet term from the EMD action. 

Consider the static spherically symmetric metrics ^, parameterized by two functions w{r) and p{r): 

ds^ = -w{r)de + + p'{r)dnl, (3) 

the scalar curvature and the Gauss-Bonnet density then read 

i? = 4j [-{^wpp' + w'p'y + 2p'[wp)' + 2] , (4) 

/:gb = 4k(V-i)]'- (5) 
p 

The corresponding ansatz for the Maxwell one-form is 

A = -f{r)dt-q.m cos 9d(p, (6) 

where /(r) is the electrostatic potential and qm is the magnetic charge. Note that the Gauss-Bonnet term breaks the 
discrete S-duality of the EMD theory without curvature corrections which is described by the transformation 

F^e-2-^*F, 0->-0, (7) 

where F — dA. It can be expected therefore that the properties of electric or magnetic black holes in this theory will 
be be essentially different. 



^ A detailed discussion of the gauge fixing can be found in jlj. In this paper we adopt the gauge gtt = l/grr- 
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A. Reduced action and equations of motion 

The one-dimensional Lagrangian associated with the ansatz is obtained by dropping the total derivative in the 



dimensionally reduced action: 

L ^ \[p'{wp)' + 1] - 2W(V' - l)'^'e'"^ - \wp^c^'^ + ipV'e^'^^ - i^e^''*. (8) 
2 2 2 z p-^ 

The corresponding equations of motion read: 

4aa[(V'-l)0'e'°l'-4aawV'0'e'''^-PP"-p'0" = 0, (9) 

4aa(Wp>'e2'''^)' - iw'V ~ ( V)' - u^P*^" + P/'^e'"^ + 4-e'°'^ = 0, (10) 
2 p-^ 

(V'^')' + 2aaK(V'-l)]'e'"* + ap'/"e2"^-a^e2°'^ = 0, (11) 

P 

(p2/'e'"^)' = 0. (12) 

Integrating once the last equation for the form field (fT2| . 

f\r) = qep-^e-^'^t (13) 



where is the electric charge parameter, one can then insert (I13p into the previous equations to obtain the six-order 
system consisting of three differential equations of the second order with the electric and magnetic charges qe , q-m 
entering as fixed parameters. 



B. Symmetries and conserved quantities 

The action ([T]) is invariant under the following two-parametric group of global transformations: 

w^we-^\ p-^pe\ (j)-^ <!)+-, f^fe-^^; r^r + v, (14) 

a 

which generate two conserved Noether currents 

drJg = 0, (15) 

9=0 

where stands for the set w,p,(j),f, and g = S, v being the transformation parameters. The conserved quantity 
corresponding to v is the Hamiltonian 

H = lip'iwp)' - 1] - 2aaw'{3wp'^ - l)0'e2'^* - ^ V^" + JpV'e''^^ + l^e'^^'t (16) 
2 2 2 2 

This quantity must vanish on shell for diffeomorphism invariant theories, so = 0. The Noether current corresponding 
to S leads to the conservation equation drJs — with the current 

Js = - ^ + 2qJ + 2a [(V' - " ^awcj)') + 2aww' pp'cj)'] e^""^ , (17) 

which is an Abelian counterpart of the integral given in [6^ . 

The integrals of motion allow one to reduce the order of the system by two leading to the forth order system. 
Moreover, for q^ ~ Q one can further reduce the order to three introducing the new variables 

w^exp(w), p^exp(^p-^^, ^ " '^^^^ 

Using them we can exclude from the system w, while w' and w" still remain. For numerical integration we will still 
use the initial six-dimensional system, applying the integrals of motion to control accuracy of the calculation. 



dgr 
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The symmetry group is enlarged on-shell on the space of the solutions of the equations of motion. It can be easily 
seen that the solution space is invariant under a three-parameteric group of global transformations which consists in 
rescaling of the electric charge 

Qe^qee^^, qrn^qm, (19) 

(leaving the magnetic charge invariant), rescaling and shift of an independent variable 

r ^re'i+^ + 1^, (20) 
and the following transformation of the field functions: 

w^we'^, p^pe', (t>^<t>+-: f^fe^-'- (21) 

a 

These do not leave the action ([1]) invariant, unless a condition 

fi = -2S (22) 

is imposed in which case we go back to the transformations (|14l) . However, dealing with the solutions, we can assume 
both parameters n and S to be independent. 



III. LOCAL SOLUTIONS WITH AdS2 x 5^ HORIZON 



We are looking for extreme black holes carrying both electric and magnetic charges in the EDGE theory for which 
the metric function w(r) has double zero at some point r — rn (the event horizon) and is non-singular for r > ru- 
Already in the EMD theory without curvature corrections the analytical dyonic solutions are known only for two 
special values of the dilaton coupling, namely a? = 1 and = 3 (for a more detailed discussion see Appendix A) 
while for generic values of a only a numerical analysis is possible. So we can hope to solve the problem of dyonic 
black holes in the EDGE theory only numerically. Meanwhile, already from a local analytical analysis of the solution 
near the event horizon we obtain important restrictions on the parameters. 



A. Horizon expansions 



Local solutions in the vicinity of the event horizon r = rjj can be constructed expanding them in terms of the 
deviation x = r — rn'- 



oo 



"(0=E""=^'' p(r•)=^Pfca;^ P(r) -e^'^^M -^Pfea;^ (23) 



k=2 k=Q k=0 



According to the assumption of extremality, the function w{r) starts with the quadratic term: vanishing of wq means 
that r = r/f is a horizon and vanishing of wi means that the horizon is degenerate. Such an expansion contains only 
one free parameter Pi for fixed values of charges: 



pI Qaa^pt 
Pi 
Aaa'^PQ 



p{r) = po + [(a' - l)qtn + M^a" ~ 2)ql + Aa^a^ - I)] x + Oix^), 



The physical value of the horizon radius po is not a free parameter here: it is fixed by the charges as follows 

2 2ge(2a + ggJ 
Po = — /,,■} ■ (25) 

Note that the dilaton coupling constant enters the expansions only through , so the space of solutions is symmetric 
under a — s- —a, (j) — > —0 (we assume a > 0). When a ^ 0, some of the expansion coefficients diverge unless Pi = 0. 
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If Pi = 0, the higher order coefficients in the dilaton expansion also vanish and the dilaton is constant. Then the 
Gauss-Bonnet terms in the action becomes a total derivative, so the theory reduces to the Einstein-Maxwell one. One 
can then show that the series expansions for w and p then combine indeed into the extremal Reissner-Nordstrom 
solution. 

Higher order expansion coefficients (see Appendix B) contain in denominators a sequence of the following combi- 
nations of the parameters 

T,(a, g„) := (a^ - a^^)qi + 2ai3a^ - 2a1)ql, + Aa'ia'iat + 2) - a^], i > 2, (26) 

where of are the integers: 

«? = l + 2 + •••+^=^^. (27) 

Therefore the expansions do not exist for the values of a and satisfying the equations 

T, = 0, z > 2, (28) 

which define the sequence of limiting curves in the parameter plane a, q„i. The only possibility to avoid the divergences 
of the expansion coefficients would be to choose the parameter Pi as the product Pi = c 11^=2 with finite c. However, 
in this case there will be a turning point outside the horizon in which dp/dr ~ 0, so that the solution does not extend 
to infinity but ends up in a singularity. 

Remarkably, the sequence of integers ((27|) coincides with that found by Poletti et al. [13] as condition of existence 
of asymptotically flat dyons with two (non-degenerate) horizons in the EMD theory without curvature corrections. In 
our case this sequence enters in the definition of the set of functions Ti(a, Qm), whose vanishing marks non-existence of 
the extremal local solution. In the limit qm ^ oo which is equivalent to a — > (decoupling of the Gauss-Bonnet term) 
the solution of the equation Ti(a, qm) = is just a = ai. Specializing the path in the a, qm plane one can in principle 
make higher order coefhcients non-singular, but the procedure is somewhat subtle (a more detailed discussion will 
follow). 

Another characteristic curve in the parameter plane a, qm is defined by vanishing of the linear term in the expansion 
of the radial function p{r): 

Ti := (a^ - l)q^ + 2a{3a^ - 2)q^ + Aa^ia"^ - 1) 0. (29) 

This signals a potential singularity of the corresponding global solution (the horizon itself is a turning point for the 
radial variable). Note that the curve a{qm) obtained as the solution of the equation Ti = reaches the value a = 1 
for qm = 0, and qm — > oo, and it has a local minimum = 4/5 for qm — \f2a. It is worth noting that the expression 
for Ti does not follow the general formula for (valid for i > 2). 

Therefore, the family of curves Tj = divide the two-dimensional parameter space of a (vertical axis) and qm 
(horizontal axis) of the global solutions for a fixed a into the disconnected regions. This situation will be described 
in detail in Sec. IV. 

The values of the integrals of motion corresponding to the series expansions (j24p are 

[d^O + qlP^^ - pV\ = 0, -h - 2 q, /o, (30) 

where Pq is given explicitly in and /o is the value of the electrostatic potential on the horizon. 

The unique free parameter Pi in the near horizon expansions will be fixed by asymptotic flatness. This means 
that the extremal dyonic solutions are completely characterized by the charges. It is convenient to absorb the Gauss- 
Bonnet coupling constant a. into the redefinition of the parameters q^ = (iel \foi^ Qm — \/oiQm and Pi Pi /a. Then 
the near horizon expansions then will read: 



"'W = ^-^7^[3(a2-l)g;^ + 6(3a2-2)g^+4(5a2-3)]x3 + 0(a;4), 



p{r) = PO + [(a' - l)qt + 2(3a' - 2)ql + 4{a^ - 1)] x + 0{x'), 



4a2po 



o2 



P(r) = ae'^^ = +Pix + Oix% (31) 
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and the relation (l25l) becomes 



pl = ^A^^. (32) 

It is important to determine the correct sign of Pi. To be able to interpret the region r > r// as an exterior of the 
black hole, one has to ensure positiveness of the derivative p' at the horizon. From the near horizon expansion of p 
one finds 

p'U=o = T^ti > 0, ti := (a2-l)g,f„ + 2(3a2-2)C + 4(a2_i). (33) 

Thus, we should take positive Pi for Ti > and negative A for Ti < 0. Introducing the sign parameter ? = we 
find therefore: 

ti 



|Ti| 



(34) 



Another useful redefinition is based on the observation that po and Pi enter the expansions in the combination 
b — |Pi|/(a^Po)- Consider now the transformations of the expansion parameters under the symmetries of the solution 
space [^T|) . It is easy to see that the full set of local solutions can be generated from one particular solution with 
P{) — 1^ b — 1^ which we will call the normalized local solution, by the symmetry transformations with 5 — — \n p^ and 
p = 21n(5po)- The normalized local solution does not contain free parameters (for fixed qm)'- 

w[r) = x2 - [3(a2 - l)gf„ + 6(3a2 - 2)ql + A{ba^ - 3)] + 0{x^), 
p{r) = l + ^-^[{a^-l)q^^ + 2{ia^-2)ql + A{a'-l)\x + 0{x^), 

= + (35) 

2(2 + g4) 

Note the presence of the sign function <; in the odd power terms. The electric charge corresponding to the normalized 
local solution is given by 

. V4 + ql. 



B. Reissner- Nordstrom limit 



One case in which the extremal Reissner-Nordstrom solution is valid is the already mentioned limit a — s- 0. However, 
there is another limit in which our local solutions make contact with the Reissner-Nordstrom solution: qm 00. This 
can be implemented by taking either 5^ — > 00 or a ^ 0. In this limit, the near horizon expansions (j3ip reduce to 



w{r) 
P{r) = Po 



^ [1 - 2Fix + 3F?x^ - 4F?a;3 ^ 0(x'^)] 
Po 



1 + Fix- 



a'-2 1 (Pig 
2a2 qi^ 



pI 



2 K-l)(2a* + 19a-'-36) 1 { Piq^„ 



24a4(a2 - 6) g; 



r,4 



Pi 



x^ + 0{x^) 



P{r) = ae^-^-^ ^ ^ + 

2 9 m 



X - Fix^ 



(a2_ 1)2(^2 +4) p^q4- 



16a4(a2 - 6) 



Pi 



+ Oix^) 



(37) 



where 



^^1 



(38) 
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This is valid for all values of a except a = Qi located on the curves of = when qm ^ oo. To ensure the 
AdS2 X structure of the horizon and the regularity of expansions we assume po to be finite and impose the 
condition Fi — const 7^ 0. Since p§ = 2qeqm, this implies Pi ^ ^Jaq^ — > 0. Then, the expansions simplify and admit 
summation to the closed expressions 

w{r) ^ —{1 + Fix)-^ , p{r)^ po{l + Fix), P{r) ^ Pq. (39) 
Po 

Now, the asymptotic flatness means -Fi = 1/po, and finally assuming rn = po we will have x = r — po which lead to 
the Reissner-Nordstrom dyonic black hole 

w{r) (1 - ^) ^ , p(r) = r, P(r) = Pq = ae""^" , (40) 

with the electric and magnetic charge parameters (je — \f^Po and = po/\/2Pq (corresponding to equal dilaton- 
rescaled charges Qe = Qm = Po/ V^) . 

Note again that there are two different ways to implement the limit q^n ^ 00: either a — s- 0, Pq ^ and the radius 
of the horizon remaining finite po = '2qeQm = '^qelm, or g,„ — )■ 00, implying the infinitely large radius of the horizon. 



C. Relation to EMD dyons 



The dyon solutions of the EMD theory without curvature corrections for a generic a exhibit one horizon and does 
not admit an extremal limit. For a discrete sequence a = Ui there are solutions with two horizons which may have 
such limits [b^ ]. This sequence is the same as found above from a different reasoning. We expect to have a relationship 
between our solution and those of the Ref. [G^I (in the extremal limit) when q„i 00 (equivalent to a ^ 0). However, 
the situation is somewhat subtle. Restarting with the set of equations of motion for a — and considering the series 
solution near the event horizon, we obtain in the lowest order: 

^0 = — , pI = 2qeqm. (41) 

qm 

This corresponds to the limiting form of the coefficients in the EDGE theory. Analyzing the higher-order equations 
for the expansions coefficients we find the following. In general pi is a free parameter (which has to be fixed by the 
asymptotic conditions) while all Wk and pk with k > 2 are completely determined by the equations order by order. 
However, when one tries to solve the equations with respect to Pk for fc > 1, an interesting bifurcation behavior is 
observed. There are two possible cases: either = 1 (then Pi is a free parameter) or Pi = (we leave aside the 
special case a = 0.) In the first case, all higher order Pk are fixed by the equations. In the second case we observe 
another bifurcation: either = 3 (and P2 is then free) or P2 — 0. Again, in the first case P{r) is fixed by P2, but 
in the second case we have a further bifurcation: either = 6 or P3 =0. The analogous bifurcations exist in any 
order. This branching procedure reproduces the value a = at i-th step. This indicates that the extremal dyonic 
black holes can exist only for this discrete sequence. 

Thus, for any i > 1, there are two independent parameters pi,Pi in the local solution {po,Pq being related to the 
charge parameters), and the expansion coefficients pj, j = 2, . . . , 2i — 1 and Pj, j ~ I, . . . ,i — 1 are all zero. Moreover, 
the expansion of z«(r) differs from the corresponding expansion of the Reissner-Nordstrom solution (a ~ 0) 

wMr) = , / =\x'^^x' + ^x'-^x' + ^x<^ + 0{x'), (42) 

(po + pixy pf, pi^ p* pf; pg 

only starting with the term W2i+2- 

In view of such a behavior, to reach the series expansions arising in the EMD theory, which could be expected to 
arise in the limit a = of the EDGE model is somewhat problematic. First, in the EDGE case there is only one free 
parameter. Pi, while there are two parameters, pi and P^, in the EMD theory. So the solutions emerging in the limit 
— *■ 00, if exist, can contact the corresponding solutions in the EMD theory only for a special value of pi. Secondly, 
the limit a — > a^, g„i ^ 00 hi the parameter space depends on the direction chosen. In particular, one gets essentially 
different results taking first a — *■ and then ^ c», or first qm 00 and then a ^ ai. And both these two do not 
seem to give the result of the EMD theory. 

We refer the reader to Appendix A for details concerning the cases of lower values of the sequence for which a 
closed form summation is possible. 
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D. Entropy and temperature 

Following the Sen's entropy function approach, the entropy of extremal dyonic black holes can be straightly calcu- 
lated 

S = 2nq,^q^:;TA^^^pl + ^^^. (43) 

Technical details are similar to those in the pure electric case treated in There are two interesting limits of the 
above expression. Firstly, the Bekenstein-Hawking entropy-area relation, S = A/A, A — Attpq, is recovered when 
a = or — > oo. Secondly, if the magnetic charge parameter is vanishing, we recover the result obtained for the 
pure electric case [ll^3|' namely, the double Bekenstein-Hawking value. Such entropy enhancement for small black 
holes was discussed in Ref. [11] ■ However, for a generic extremal dyonic solution, the black hole entropy can not be 
completely expressed in terms of its horizon area. 

It is worth noting that although the value of the entropy can be calculated using only the local solution valid in 
the vicinity of the AdS2 x event horizon, its interpretation as the entropy of a black hole presumes an existence 
of the global solutions extending to infinity. We will see later on that this imposes certain restrictions on the values 
of the magnetic charge and the dilaton coupling constant a. Purely local analysis is therefore insufhcient for drawing 
conclusions about the correspondence between the string and geometric values of the entropy. 

The temperature of the extremal EDGB black hole is zero, as for the extremal solution without the Gauss-Bonnet 
term: 



1 / , dy^t 



27r V^^ dr 



^{r-rH)\r=r^^O. (44) 



IV. ASYMPTOTIC BEHAVIOR 



Another region where local solutions can be constructed analytically is the asymptotic zone r —> oo. Similarly 
to the case on the uncorrected EMD theory, we find that two type black holes can exist: usual asymptotically flat 
solutions, and black holes on the linear dilaton background [65]. In the latter case the dilaton diverges at infinity 
(linearly for a proper choice of the radial coordinate), but the ADM mass of the black hole itself is finite, so that the 
solution can be interpreted as the black hole on the linear dilaton background. 



A. Asymptotically flat solutions 

Looking for asymptotically flat global solutions we have to ensure w —^ 1, p/r 1, cj) ^ (j)oo (constant) as r — > oo. 
The subleading terms should be expandable in the power series of 1/r. The asymptotic solution with these properties 
contains five parameters: the ADM mass M, the electric and magnetic charges Qe, Qm (rescaled), the dilaton charge 
D, and the asymptotic value of the dilaton (j)^: 

n^ir) = l-^ + ^^^i±^ + 0(.-), 

2r 3r^ 
D 2DM - aaQl + a'^aQ 

where 



^qee"^"^-, Q^^g^e-^"^-. (46) 
The dilaton charge can be also read off from an asymptotic expansion of the dilaton exponential: 

The values of two integrals of motion in terms of the asymptotic parameters are: 

^ = ^(^oop'c^-l), Js = 2q,J^-M~-. (48) 
2 a 
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The constant /oo in the electric potential is usually fixed to zero. 

Behavior of the global solution which starts with the normalized local solution ([55]) at the horizon depends only on 
the dilaton coupling constant a and magnetic charge parameter — \/aqm. Its existence for all a is not guaranteed 
a priori. But, in some intervals of a whose boundaries depend on q„i, as can be shown numerically, there exist 
solutions varying smoothly with increasing x such that the function w and the derivative p' stabilize at infinity on 
some constant values Woc 7^ 1, p'oo 1- Then, using the symmetries (|^ [^T|) of the solution space, one can rescale the 
global solution obtained to achieve the desired unit values for these parameters. As we have argued, two parameters 
fi,6 effectively replace the parameters po,Pi of the (non- normalized) local solution or (|3ip . So one could expect 
that rescaling of the solution ensuring Woc — 1, p'oo — 1 would fix both quantities po, Pi on the horizon. But from the 
Hamiltonian constraint equation H = with H given by the Eq. ([^5]) it is clear that one must have Woop'^ = 1 for 
any asymptotically flat solution. Therefore it is enough to perform one but not two independent rescalings in order 
to achieve Wao — 1, p'oo = 1- Indeed, under the transformation (j201 12ip . the relevant functions and parameters are 
transformed as follows 

w^we", p'^p'c^/^ wp'^^wp'^: po ^ poe\ Pi^PiC^-f'^. (49) 

Since the choice of /i, 5 is equivalent to the choice of poj Pii an invariance of the product wp'^ means that the solution 
starting on the horizon with any po,Pi will reach at infinity the values Woo, p'oo satisfying Wq^p'^o — 1- Therefore, 
taking p — —\vl Woo , we will achieve simultaneously Woc — 1 and p'^ — 1 . This means that asymptotically flat 
solutions still form a two-parameter family, two parameters being the electric charge q^ and the magnetic charge q„i. 
Five asymptotic parameters M, D, (j) 

00 ■) Qe-) Qm Q-rG functions of ci^y^ which can be found numerically. 



B. Black holes on the linear dilaton background 



There is another type of physically interesting black hole solutions which asymptotically approach the linear dilaton 
background (LDB) [63]. For such black holes, some metric functions div erge asymptotically, so to be able to recognize 
them numerically we have to pass to some conformally rescaled metric [69| 

rf4ual = e-2'^W-*-)d.^ (50) 
which has the following explicit form in our case 

In the dual frame the asymptotic metric for LDB solutions could be either x S"^ (only for a special value of the 

dilation coupling, a? = 1 for the 4-dimcnsional theory) or 74^5*2 x S*^ [H, [tO]. Therefore, if we rewrite the metric in 
the following form 

ds^^^j ^ -pdt^ + ^du^ + p-^p^dn^, (52) 

where du = dr/ P, then the metric functions should have the following asymptotic limit (determining the radii of 
the S^) 

P-'p^ -Ro, ^-Rlu\ (53) 

where i?o and i?o denote the curvature radii of the 5^ and AdS2 respectively. The differentiation of the second 
equation gives the relation 

which can be used in the numerical procedure. 



V. NUMERICAL ANALYSIS 



In this section we present numerical results for the global extremal dyonic black hole solutions. We extend the local 
solutions constructed via series expansions near the horizon to the asymptotic region by numerical integration. As 
expected, the free parameter Pi turns out to be fixed by the asymptotic conditions, either of asymptotic flatness or 
the linear dilaton asymptotic (LDB). 
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A. Asymptotically flat dyons 



For the pure electric extremal black holes pj, the global solutions were found to exits only in a limited range of 
the dilaton couplings less than a critical value acr- By turning on the magnetic charge, the range of a for domain 
of existence can be extended from the interval < a < Ocr to an infinite sequence of disconnected intervals, located 
between the limiting curves = in the two-dimensional parameter plane a, qm- More precisely, the regular solutions 
exist for a satisfying 

a~{qyn) < a < af{qm), (55) 

where a~ is located above the curve Ti_i = and dl below the curve = 0. The domains of existence obtained 
numerically are shown in Fig. 1. The most surprising feature of this plot is that for q^n oo, which is equivalent 
to turning off the Gauss-Bonnet term, the domain of existence of dyon solution does not reduce to discrete values 
of the dilaton coupling constant a = a^. As we have seen, already the series solutions in the near-horizon region 
are essentially different in the EDGB model and in the curvature uncorrected EMD theory: the former being one- 
parametric, while the latter — two-parametric. Transition between these two series solutions is unclear, the limiting 
form of the series in the EDGB model depends on the direction in the a, q^ plane in which the limit qm — *■ cxd is 
taken. Numerical solutions were checked to exist up to some large values of qm, that is, with small Gauss-Bonnet 
coupling. They still exist for continuously varied a, not for discrete values as in the pure EMD theory. Therefore the 
curvature corrected theory gives qualitatively different predictions for dyons even when the Gauss-Bonnet coupling is 
small. 

Our numerical results reveal the following behavior of dyon solutions for large q^, (or small Gauss-Bonnet coupling). 
If we fix qm and then take the dilaton coupling a close to the critical values af, the parameter Pi (which is fixed by 
the asymptotic flatness conditions Woo ~> liP'oo ~* 1) going to diverge, and the solution becomes ill-defined. 

However, if we fix the value a ^ of the dilaton coupling and then move g„j to infinity, the parameter Pi approaches 
zero. In this case, the limiting solution qm oo corresponds to the dyonic Reissner- Nordstrom extremal black hole 
(with the frozen dilaton). Such solutions, however, do not exist in the pure EMD theory. Moreover, the subsequent 
limit — > of always gives the dyonic Reissner-Nordstrom solution. Therefore, the extremal dyons in the pure EMD 
theory look like a set of discrete points in the two-dimensional continuum which is difficult to resolve. Nevertheless, 
the limiting value of the radius of the horizon is the same both in the EDGB and the EMD theories (see also Appendix 
A). ^ 

Using the symmetry of the solution space under the ^-transformation (jl9[) , one can generate the sequence of solutions 
with different electric charges qe and correspondingly with different masses, dilaton charges and the asymptotic values 
of the dilaton (j)ao ■ Since variation of the electric charge is essentially equivalent to variation of the unique parameter 
po (for fixed qm) in the horizon expansion, it is clear, that using (5-transformation we will generate all extremal 
solutions. Under this transformation the mass and the dilaton charge scale as e'', while the electric charge and the 
dilaton exponent e^"*^" scale as e^^ . Therefore the ratios 

kM = , ko = — , few, = , (56) 

qe qe qe 

depend only on a. Their numerical plots are presented in Fig. 2. 

With growing magnetic charge, two end points of the allowed interval of a tend to the boundary curves. When the 
dilaton coupling approaches the critical values (the end points of of each segment of the regular solution with fixed 
qm), the global physical quantities such as the mass M, the dilaton charge D, the electric charge Qe — qe 6^°"^°° and 
the magnetic charge Qm = qm '^'^ diverge and the only free parameter Pi also goes to infinity (in order to ensure 
an asymptotic flatness). However, the ratio of the following physical quantities 

has a simple limit ksps 1 (Fig- 3). Remarkably, as a ^ a^, the associated global parameters tend to satisfy the 
following relation 

aHP + D' = — {Ql + Ql\ (58) 

which coincides with the BPS condition for charged black holes in the EMD theory without curvature corrections. 
One can also see that in this limit the discrete S-duality of the EMD theory is restored. This feature is similar to 
that in another stringy generalization of the EMD theory in which the Maxwell action is replaced by the Born-Infeld 
action but no Gauss-Bonnet term is introduced 17111. 
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FIG. 1: The domains of existence of EDGB dyons in the two dimensional parameter plane a, qm for qm < 6. The family of 
limiting curves corresponds to solutions of the equations Ti = implying singularities of higher order coefficients in the horizon 
series. Black circles correspond to asymptotically flat solutions, red squares — to black holes on the linear dilaton background. 




FIG. 2: Asymptotic parameters AI, D, 0oo in terms of the horizon parameters qe, qm as given by the ratios /cm — ko ~ — 

2a^£> 

and — depending on qm ■ These ratios are shown as functions of the dilaton coupling constant a for some values of the 

magnetic charge: Qm = 1 (black, thick), = 2 (gray, dashed) and qm ~ 4 (red). 



B. Magnetic black holes on the LDB background 



The domain of existence of black holes on the linear dilaton background satisfying the condition ([55]) is located to 
the left from the condensing family of limiting curves Tj = on the parameter plane a, qm in Fig. 1. The boundary 
curve is given by Too = or, explicitly, 

qi+Aaql-Aa^{a^-l)^0. (59) 

Domains of existence of all asymptotic flat solutions are located to the right from this curve while those with the LDB 
asymptotics — to the left. Asymptotically LDB dyons do not exist for dilaton couplings less than one. Moreover, for 
larger magnetic charges, the asymptotically LDB solutions require larger value of the dilaton coupling. 

From an analysis of similar solutions in the EMD theory [t^I it follows that one of the charges of the dyonic 
configuration defines the strength of the background electric or magnetic fields in the LDB, while another is associated 
with the charge of a black hole. So one can have electric black holes on the magnetic LDB or magnetic black holes 
on the electric LDB. The asymptotic form of the metric is 

p^(eor)5, (60) 

where is the scaling parameter associated with the field strength of the background, and the dilaton for the electric 
black hole on a magnetic background is 



(61) 
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FIG. 3: The ratios ksps 
(gray, dashed) and - 
domains. 




while for the magnetic black hole on an electric background 




(62) 



The solutions presented in this paper are consistent with the dilaton field ([5^ and the scaling parameters /i = — 2(5 = 
— In^o- Thus, our solutions can be interpreted as magnetically charged black holes on the electric LDB (recall that 
the discrete S-duality is broken in the EDGB theory). In the case of the vanishing magnetic charge they are physically 
expected to reduce to the pure electric LDB without a black hole. Technically, however, our solutions can not have 
such a limit since we have assumed the existence of the horizon a priori. 



In this paper, we have shown that the EDBG four-dimensional gravity admits extremal dyonic black hole solutions 
with the horizon of the AdS2 x S'^ type. Somewhat surprisingly, adding the Gauss-Bonnet term to the Einstcin- 
Maxwell-dilaton theory leads to an enhancement of the domain of parameters for which the global solutions exist. 
Namely, the asymptotically flat dyon solutions in the EMD theory exist only for a discrete sequence of the dilaton 
coupling constant values, while in the model with the Gauss-Bonnet term the continuously varying parameters are 
allowed. An effective parameter space is a two-dimensional plane which is split into the sequence of regions separated 
by the limiting curves, marking singularities of the coefficients of the local power series solutions. These curves are 
related to the above discrete sequence of parameters of the EMD dyons, and when approaching them, the solutions 
of the EDGB model exhibit saturation features similar to the BPS conditions of the EMD theory. 

The relationship between the extremal EDGB dyons and those in the EMD theory is non-trivial. Dyon solutions 
of the EDGB theory exist only with non-zero electric charge, and for large values of the magnetic charge they tend 
to the Reissner-Nordstrom solution with a frozen dilaton and not to the discrete family of extremal dyons in the 
EMD theory as could be expected. Therefore, the pure EMD theory without curvature corrections predicts different 
black hole solutions than the corresponding curvature-corrected theory in the limit of vanishing curvature coupling. 
Remarkably, the latter limiting theory leads the same value of the horizon radius, as the pure EMD theory. 

The entropy of the extremal EDGB dyons interpolates between the Bekenstein-Hawking value in the limit of the 
large magnetic charges (equivalent to the vanishing curvature coupling) and the doubled Bekenstein-Hawking value 
in the limit of purely electric solutions. The entropy can be calculated using only local analytical solutions valid in 
the vicinity of the event horizon. The expression obtained, however, does not bear any sign of the bounds on the 
parameters for which the global black hole solutions exist. Such bounds can be revealed once we try to extend the 
local solutions to infinity as asymptotically fiat ones. This lesson is worth to be kept in mind in the discussion of 
the entropy of the curvature corrected black holes in string theory. Purely local analysis of classical solutions is still 
insufficient to fully understand the entropy of black holes. 

We also found the second family of the EDGB dyons which asymptotically approach the linear dilaton background. 
For them, the metric is not asymptotically flat, but physical value of the black hole mass and other parameters can 
be extracted by subtracting the background values. In string theory such black holes correspond to thermalized 



VI. DISCUSSION 
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states of the quantum theory through the QFT/Domain wall correspondence. Our solutions may be interpreted as 
magnetically charged extremal black holes on the electric linear dilaton background. 
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APPENDIX A: EXTREMAL DYONS IN 4D EINSTEIN-MAXWELL-DILATON THEORY 

The action and the equations of motion of the four-dimcnsional Einstein-Maxwell-dilaton (EMD) theory with 
arbitrary coupling a are obtained setting a = in the Eqs. ([51 llip . Consider again the series expansions in the form 
(I24p around the horizon r = ru in terms of the deviation x r — rn ■ In the lowest order we obtain: 

Pq^ —, pI = ^qeqm, (Al) 

q,n 

which coincides with the a limit of the corresponding EDGE relations. However, calculation of the higher 
order coefficients shows that one deals with the local solution containing two free parameters (with fixed electric and 
magnetic charges), which can be conveniently chosen as pi and Pi. Recall, that the corresponding local solution 
in the EDGE model constructed in Sec. 3 contained only one free parameter Pi. Another new feature consists in 
bifurcations encountered in the calculation of the higher expansion coefhcients as was described in Sec. 3. These 
bifurcations give rise to a sequence of discrete values of the dilaton coupling constant af = i(i + l)/2 for which 
only one obtains the non-trivial expansions. For each i the series expansion for the dilaton function starts with the 
term of the i-th order. Pi, with all lower coefficients being zero, Pj<i — (except for Pq). This discrete sequence 
was also observed in the Ref. [G^I where it was shown that the asymptotically flat solutions exhibit two horizons if 
a = 0, 1, a/S, ■ • ■ , \/ n{n + l)/2, and one horizon otherwise. Since the solutions with the second order zero of the 
metric function gtt at the horizon can be obtained only starting with the non-extremal solutions with two horizons, 
this is consistent with the special values we find here for genuinely extremal solutions. Note that only in the first two 
cases a = 1 and a = \/3 the exact dyonic solutions have been obtained analytically. In what follows we show how this 
can be done via direct summation of the series solutions. 



The series solution reads 



pI pI V Po Po / V Po pI J V Po Po Po° J 

6^^20p3p,V.^6p^X 

Po Po Po / 

^llni ^2 , Pl^l^m^S Pl^llrn , ^1 1m \ 4 , Pl^llm , ^Pl^l1rn\^b 



„4p2 4 firt2p4 8 p6„12 
Pl^l Im , "Pl^l 1m , A 1m 



2pl 8p9 16pii 



O(x^), (A3) 



Pi Pl1m\ ^2 , D /^Pl Pl1m\ „3 u Pllm \ 4 



P(r) = Po + Pix-Pi[^- 1^ a;^ + ^li _ llllli - Pi I ^ - x^ + 0{x^). (A4) 

\Po Po / VPo Po / Vpo Po 
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The summation is performed by collecting terms of different order in pi inside the series in terms of P^q^- For 
example, for p(r) one can rewrite the series as 



= (po + pi '■'■) 1 



p2 4 



1 Pi , Pi 2 

1 X H — ttx 

Po Po 



p4 8 



Po 



^1 1: 



X 



+ 



p4 8 



1 



2P(1 (Po + Pixy ' (po + Pia;)^ 

= ^(po + Pia;)2-P^^a;2. 

A similar pattern can be found in the series expansion for w{r), the result being as simple as 



(A5) 



w(r) 



p2(r) ■ 



For P(r), introducing a new parameter Pi defined via Pi = PiPo one gets just a geometric recurrence: 



P(r) = Po + Pia; 
= Po + 



Pi Piql 



Po 



p§ 



Pi PioL 



Po 



P§ 



pgPia: 



= Pr 



(po + pix) + Pig^x 
(po + Pi.^;) - PiqliX " (po + pix) - Piql^x 

Thus, we have obtained an exact solution in a closed form. 



Pi Piqi 



Po 



Pi 



x' + .- 



(A6) 

(A7) 
(A8) 



2. a^ = 3 



This case corresponds to the Kaluza-Klein theory. The series expansions are 



^(r) - la:'-^x^^+'^x^-^x'+(^-^ + ^-^\x<^-(^-^ + ^P^ 

^ ^ " Pi Pi Pt Pi \pi 9pe r \pI ^ H 

32 „4 



P^ 
32 „4 



Po 9p^ 

,2 p2„4 



Pi 



x^ + 0(x8). 



P2„4 ^ p2„4 0o2p2„4 

p(r) = po + pra.-4%x4 + ^i|^x^-^^^x« + 0(a.^ 
ypo dpg cipo 

P(0 = Po + P..^-^.3^^(^d±3p?),4^0(.^). 

po Po 

Rearranging an expansion for p(r) in the same way we obtain the exact sum 



The function w is given again by w(r 



p(r) = 



,,.2 1^2 



(po + pix)"^ - -PiqtnX^ 



1/4 



(A9) 
(AlO) 
(All) 

(A12) 



/ p {r). However, it is more difficult to sum up the series expansion for the 



dilaton function P(r). An easier way to find P(r) is to solve the equations of motion directly using the above results 
for p(r) and w{r). This gives 



P(r) = Po 



(po + pix)2 + |P2g2^a;2T^/' 



(po + pia;)2 - gPjg^x 



(A13) 



Again, we have obtained an exact solution from the near horizon expansion. 
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3. = 6 and beyond 

One can try to repeat the same procedure for higher at. In particular, for = 6 (i = 3) one has the following 
series expansions: 

_ J_ 3,M 4_fpi 5,^1 7 

Po Po Po Po Po Po 



pI lopg; 

zu^Q 



P(r) = P, + P,x'- ^ a;^ + ^ + 0{x% (AM) 

Po Pa 

Note, that the function w{r) differs from the Reissner-Nordstrom w(r) only in terms of the order and higher 
(for higher i from the terms of the order 2i + 2). For the dilaton function an expansion starts from the cubic term 
(generally, from the i-th terms). But we were not able to find a closed form for these expansions as a whole. 



APPENDIX B: HIGHER-ORDER TERMS IN THE LOCAL SOLUTION NEAR THE HORIZON 

Here we list some higher order coefficients in the near horizon expansions for extremal black holes in the EDGB 
theory considered in Sec. 3. The parameters are po which is determined by the electric and magnetic charges, and Pi 
(free) which must be fixed by asymptotic conditions: 



Po = 

P2 = 

Ps = 



Pi 



P? 



-ia'-3){a'-l)ql 



4ap2a2[(a2 _ 3)^,4^ + Qa{a^ _ 2)^2^ + 4a2(5a2 _ 3)] 

+4a(5a^ - 6)q^ + 4a^(5a^ + 3a^ - 18)q^ - 32Q!^(2a^ + 3)g^ + 16a*(a^ - ba^ - 3) 

P[ 

48a2a4p4[(a2 _ 3)^4^ + 6a(o2 - 2)5^ + 4a2(5a2 - 3)][(a2 - 6X + 6a(o2 - 4)g^ + 4a2(8a2 - 6)] 

-3(a2 - l)2(a2 - 3)(a2 + 4)^^^ _ q^^^2 _ i)^24a^ - 91a^ - 9a'^ + 96)ql^ 

_/i^2^c.i _ i828aS + 24440^^ + 45a2 - 1008)g,\2 _ 8a^(483a^ - SOOGa^ + 6283a^ - 2145a2 - "-^ 

5625a2 - 2520)?^, + 32a-"^(165a^ - 1762, " 



- 1828aS 



-16: 



-64a^(1156a^ 



- 24440^^ ^ 

- 6647a'' 



1361a*^ 



1922a'' 



- 4785a2 



)<lni + 

- 1008)^^ 



(Bl) 



(B2) 



128. 



+512a^(23a* + 55a^ 



- 165a'' 



135a2 + 18) 



iOtia" - Lib'Za" - 1844a'* H 
!Q;^(369a^ - 49a'5 + 1573a'' 



- 6975a2 
+ 1755a2 



2016)9„ 



,10 



2016)g^ 
F 288)9^ 

(B3) 
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Pi 

P2 
P3 



Pi 



[a^ - l)qt, + 2a(3a2 - 2)g2^ + Ao?{a^ - 1) 

Pi 



2a2po[(a2 - 3)9^ + 6a(a2 - 2)q^^ + Aa'^(ba? - 3)] 



(a'-3)(a^-2)g« 



-6a(5a 



4a2(5a'' + ISa^ - 36)q^ - 8a^(5a2 - 24)qf„ - lQa^{a^ - 6) 

P3 



n 

24aoVo[(a^ - 3)g^ + 6a(a2 - 2)g^ + 4a2(5a2 - 3)][(a: 



(a^-l 



_)(a2-3)(2a4 



19a2 - 



36)9if 



- 6)g4„ + 6a(a2 
- 206a^ 



2a(66a^ - 259a^ - 



8a- 



-6)] 



12930^ 



-4a2(304a* - 1613a6 + 1085a'' + 3327a2 - 3024)gi2 8^3(4860^ - 2716^ , 

-16a''(438a*^ - 687a6 + 2704a'' + IGOSa^ - 7560)g*, - 32a^(24a* - 2273a^ + 2200a'' + 1689a' + b048)g„ 
-64a^(1098a® - 2261a^ + 4467a* + 2427a2 + 3024)g^ - 128a^(342a* - 758a^ + 2409a* + 1173a2 

-aS(9728a^ - 215040^ + 110592a'' + 53760a2 + 27648)1 ; 



r,12 



la'' 



(B4) 



(B5) 



6048)g„, 



,10 



4a2( 
864)gi* 

3599a* + 4065a2 - 

^ 1689a^ + 6048)gJ|, 

864)g^ 

(B6) 



W2 



W3 



9 ' 

Po 



Pi 



3(a2 - l)q^ + 6a(3a2 - 2)g^ + 4a2(5a2 - 3) 



P2 



48a2a4^6[(„2 „ 3)^4^ _^ 5^,(^2 _ 2)^2^ + 4^2 (5a2 _ 3)] 



9(a2-3)(a2- 1)25,1 



+ 18a(a2 - l)(9a'' - 31a2 



18)gi^? + 12a2(i06a^ + 411a2 - 396a* - 135)g^ 



+8a3(741a^ - 2150a* + 1782a2 - 540)g^, + 16a*(1010aS - 2220a* + 1413a2 - 405)?^ 
+96a^(223a^ - 404a* + 195a2 - 54)g2 + 64a^(173a^ - 269a* + 99a2 - 27) . 



(B7) 
(B8) 



(B9) 



An inspection of these and higher-order coefficients shows the systematic appearance of the combinations in 
denominators, as described in the Sec. 3. 
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